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GENE MODULATION BY RB2/i>130 EXPRESSION 



The retinoblastoma gene family consisting of RBljpl05 y 
pl07, and RB2/pl30 cooperate to regulate cell-cycle 
progression through the Gl phase of the cell cycle. 
Previous data demonstrated an independent role for the 
reduction or loss of pRb2/pl30 expression in the 
formation and/or progression of lung carcinoma. Rb2/ 
pl30 is mutated in a human cell line of lung small cell 
carcinoma as well as in primary lung tumors. To identify 
potential pRb2/pl30 target genes in an unbiased manner, 
we have utilized an adenovirus-mediated expression 
system of pRb2/p!30 in a non-small lung cancer cell line 
to identify specific genes that are regulated by pRb2/pl30. 
Using oligonucleotide arrays, a number of Rb2/pl30 
downregulatcd genes were identified and their regulation 
was confirmed by semiquantitative reverse transcription- 
polymerase chain reaction (RT-PCR) and Western blot 
analysis. As a result, 40 genes showed greater than 2.0- 
fold "modification in their expression level after the RB2f 
pi 30 viral transduction. In conclusion, coupling adenoviral 
overexpression with microarray and semiquantitative RT- 
PCR analyses proved to be a versatile strategy for 
identifying pRb2/pl30 target genes and for better' under- 
standing the expression profiles of these genes. Our results 
may also contribute to identifying novel therapeutic 
biomarkers in lung carcinoma. 

Oncogene (2003) 0, 000-000. doi:10.1038/sj.onc. 1206866 
Keywords: lung cancer; Rb2/pl30; microarray analysis 



Introduction 

Lung cancer is the leading cause of cancer death 
worldwide and one of the most common malignancies 
diagnosed in the United States. The American Cancer 
Society estimated that 169 500 new cases of lung cancer 
would be diagnosed in 2001 and that 157400 people 



would die of the disease that year. Data from the 
American Cancer Society for the year 2002 indicate that 
169400 new cases of lung cancer are estimated to occur 
and 154900 patients are expected to die because oflung 
cancer (Greenlee et a/., 2001). Additionally, lung cancer 
is usually diagnosed at an incurable stage. 

The high mortality rate for lung cancer probably 
results from the absence of effective therapies as well as 
standard diagnostic procedures of early tumoral stages 
compared with colon, breast, and prostate cancers 
(Wiest et 1997). The majority of bronchogenic 
carcinomas can be classified into four histological types: 
small cell lung carcinomas, adenocarcinomas, squamous 
cell lung carcinomas, and large cell carcinomas. Small 
cell lung carcinomas are a separate entity, whereas the 
behavior of the other three histological subtypes is 
similar, for this reason these are grouped within the non- 
small cell lung cancer (NSCLC) type. NSCLC accounts 
for nearly 80% of lung malignant tumors and it is 
associated with a poor prognosis. 

Lung cancer is the result of molecular changes in the 
cell, resulting in the deregulation of pathways that 
control normal cellular growth, differentiation, and 
apoptosis. In this scenario, proto-oncogenes and tumor 
suppressor genes are found to be mutated or have 
abnormal expression patterns in this disease. Although 
much is known about the natural history, predisposing 
factors, and outcome of NSCLC, the understanding of 
this disease is still unclear. Even though many molecular 
changes associated with NSCLC have been reported 
(Hibi et ai, 1998; Forgacs et al. 9 2001), a fully 
understood mechanism associated with this type of 
cancer has not been described yet. 

The retinoblastoma gene family consisting of RBI, 
pI07, and RB2/p/30 cooperate to regulate cell-cycle 
progression through the G I phase of the cell cycle. 
Products of RBl y pI07, and RB2jp!30 are characterized 
by a peculiar steric conformation, the 'pocket region 1 , 
which is responsible for most of the functional interac- 
tions that characterize the activity of these proteins in 
cell-cycle homeostasis (Paggi el cr/„ 1996). Rb family 
members are nuclear proteins that are regulated in a cell- 
cycle-dependent manner by phosphorylation, exhibit 
growth suppressive properties in a cell-type-dependent 



manner, are implicated in various forms of differentia- 
tion, and are critical targets for inactivution by 
transforming oncoproteins of DNA tumor viruses. 

During the past 6 years, the involvement of pRb2/ 
pi 30 in~ lung cancer has been studied. In fact, 
immunohistochemical analysis of the expression pat- 
terns of the Rb family members (pRb/pi05. pl07. and 
pRb2/pl 30) in 235 specimens of lung cancer ( Baldi et aL 
1996) and the expression pattern of pRb2/pl 30 in 1 555 
specimens of human lung cancer showed an inverse 
correlation between the histological grading of the 
tumors, the development of metastasis, and the level 
of expression of pRb2/pl30 (Baldi et aL 1997). A 
statistically significant inverse relationship between the 
histological grading and the expression of pRb/pl05. 
pi 07, and pRb2/p!30 was found in fine needle aspira- 
tion biopsies of squamous cell carcinoma patients 
(Mtnimo et aL 1999). 

Point mutation in the RB2jpl30 gene determining 
protein loss in a cell line of human small cell lung 
carcinoma (Helin et aL. 1997) and identification of point 
mutations in primary lung cancer (Claudio et <//., 2000) 
were found, indicating that pRb2/pl30 may play a role 
in the pathogenesis and progression of certain lung 
cancers. 

We have utilized an adenovirus-mediated expression 
of RB2/pI30 in the H23 NSCLC cell line to identify 
specific genes that are regulated by pRb2/pl30. Using 
oligonucleotide arrays, a number of downregulated 
genes by pRb2/pl30 were identified and their modula- 
tion was confirmed by semiquantitative reverse tran- 
scription-polymerase chain reaction (RT-PCR) and 
Western blot analysis. 

Materials and methods 

Cell lines 

The human lung adenocarcinoma cell line H23 has been 
described previously (Claudio et aL 2000). The packaging cell 
line 293 (primary embryonal human kidney cells) transformed 
bv sheared human adenovirus type 5 has also been previously 
described (Claudio ei aL 1999). H23 cells were maintained in 
DM EM supplemented with 10% fetal bovine serum. 2mM L- 
gluiamine. 293 cells were maintained in DM EM supplemented 
with 10% heat-inactivated fetal bovine serum. 2m.M L- 
glutamine. 

Adenoviruses 

Adenoviruses were generated by subcloning the full-length 
ORF of the RB2' t plW gene into the pAd.CM V-Linkl vector 
to form the Ad.CMV-RB2ipUQ virus, as described previously 
(Davis et a/.. 1998: Claudio et aL 1999). The pAd.CMV-Linkl 
vector alone (to produce the Ad-CMV virus) wis used as a 
negative control to assay the effects of viral infection alone 
without delivering a transgene. Adenoviruses were expanded, 
purified, and uttered as previously described (Claudio et aL 
1999). 

Flow cytometry analysis 

Flow cytometry analysis (FACS) of exponentially growing 
H23 cells or H23 cells transduced with Ad-CMV or Ad-CMV- 



Rb2/pl30 were carried out as previously described (Claudio 
et aL 1996). Briefly. 5 x 10* cells were seeded and 24 h after the 
cells were transduced with 50 multiplicity of infection (MOl) of 
adenoviruses. At 48 h after transduction, cells were collected 
and analysed using a Coulter Flow cytometcr. 

Microarray analysis 

Before submission of RNA samples for analysis, protein 
extracts prepared from replicate plates of the corresponding 
cell culture were analysed for expected enhanced expression of 
pRb2/pl30 using Western blots. Oligonucleotide-based micro- 
arrays were purchased from Mergen (Mergen Ltd, San 
Leandro. CA. USA). The ExpressChip H05000 DNA micro- 
array system'was used for this study. This array contains more 
than 3200 genes that are involved in a variety of different 
processes. DNase-trcated total RNA (20/<g) from H23 
(parental cells). H23 cells transduced with Ad-CMV. or Ad- 
CMV-RB2/pl30 cell lines 48h after transduction were 
extracted using TRIzol (Life Technologies. Inc.. Grand Island, 
NY, USA) according to the manufacturer's protocol. RNA 
integrity was verified for lack of degradation by formaldehyde 
eel electrophoresis. The bio tin-labeled cRNA probes prepara- 
tion, hybridization, and array scanning were performed using 
Mergen Labeling/ Hybridization/ Detection Service. Data ac- 
quisition and data analysis were performed using Imagene 
software (Biodiscovery Inc., Marina del Rey; CA. USA) and 
Mergen's ExpressDzta.™ software (Mergen Ltd, San Leandro, 
CA. USA). Briefly, data were processed for local background 
correction and normalization. The raw spot for each gene was 
calculated as the mean signal of the spot values minus that of 
the local background. The / ma , value was set to 65000. after 
local background removal. A normalization coefficient (N) 
was applied to either the control population or to sample spot 
raw values to compensate for slide-to-slide and probe-to-probe 
variations. The normalization coefficient was applied only if 
(Raw spot/ N)< / max , otherwise values were set to 1^. 
Normalized values £0 were excluded from the analysis. Genes 
regulated by adenovirus transduction (Ad-CMV) with respect 
to'the parental cell line were removed from the analysis. Spots 
with mean intensities > 45 000 were excluded for the ratio 
analysis. The expression ratios calculated with corrected values 
less than mean of the local background on both channels were 
not used. Expression . ratios of the analysed genes were 
calculated comparing genes* expression values of H23 cells 
transduced with RB2/pl30 with those.of parental H23 or H23 
cells transduced with Ad-CMV. A twofold or higher levels of 
target genes* expression ratio was considered significant, in 
accordance with most of the literature. 

Northern blot analysis 

H23 cells were grown to 60% confiuency. then infected with 
50 MOl of adenoviruses carrying the RB2/pl30 gene or with 
the control Ad-CMV. After I4h, the medium was changed, 
and cells were harvested at 48 h from the transduction. 
DNAse-treated. total RNA from H23. H23-Ad-CMV, and 
H23-Ad-CMV-Rb2/pl30 transduced cells were extracted using 
TRIzol (Life Technologies. Inc.. Grand Island. NY. USA) 
according to the manufacturer's protocol. Northern blot 
analvsis was performed as previously described (Claudio 
ei aL 1994). 

Semiquantitative RT-PCR 

RT-PCR was used to analyse target gene expression in the 
presenr study. A 2/<g aliquot of DNAse-treated total RNA 
from each sample was reverse transcribed for single -stranded 
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cDNAs using tM-MLV reverse transcriptase (Invitrogort. 
Carlsbad. CA. USA) according to the manufacturer's proto- 
col. The same cDNA product obtained from each sample was 
used for subsequent PCR amplification with the primer sets 
prepared for the target gene and jUtciin [Act-li)/HPRT 
housekeeping genes. The amplification of the fi-actin and 
HPRT genes were used as double internal control. The ratio 
between the samples and each housekeeping gene was 
calculated to normalize for initial variations in sample 
concentration and as a control for reaction efficiency. Primer 
sequences were designed using the software Primer 3 (devel- 
oped by Steve Rozen. Helen J Skaletsky) available on-line at 
http://www-genome.wi.mit.edu. Primer sequences can be 
provided upon request. PCR reaction conditions were 
individually optimized for each gene product studied and the 
number of PCR cycles was set up to be within the linear range 
of product amplification. 

In each experiment, possible DNA contamination was 
determined by a control reaction in which reverse transcriptase 
was omitted from the reaction mixture. PCR products were 
loaded onto ethidium bromide stained 1.5% agarose gels. 
Densitometry analyses of the PCR products were performed 
using an Alpha Imager system (Alpha Innotech Corporation, 
San Leandro, CA, USA) and the ImageJ vl.29 software 
(developed by Wayne Rasband) available on-line at hup:// 
rsb.info.nih.gov/ij/. All PCR products were purified using 
QIAquick PCR purification kit (Qiagen. Santa Clarita, CA. 
USA) and their identities verified by automated DNA forward 
and reverse sequencing using a dideoxy terminator reaction 
chemistry for sequence analysis on the Applied biosystem 
Model 373A DNA sequencer. 

Western blot analysis and antibodies 

Western blot analysis of exponentially growing H23 cells or of 
H23 cells transduced with Ad-CMV or Ad-CMV-Rb2/pl30 
were carried out as previously described (Claudio ei aL 1999). 
Extracts were normalized for protein concent by Bradford 
analysis ( Bio- Rad Laboratories. Inc.. Melville. NY. USA) and 
Coommasie blue sel staining. 

Primary anti-B-MYB. E2F-1. KPNA2. MK.K3. NIK. 
PCNA. PLK, RAFI (Santa Cruz Biotechnology. Inc.. 
Santa Cruz, CA ? USA). ahti-MAGE-A (Upstate. Lake 
Placid, NY. USA), and anti-HSP70 (Oncogene Science. 
Cambridge. MA. USA) were used following the manufac- 
turer's instructions. 



Results 

Effects of RB2/pJ30 adenoviral transduction on the H23. 
lung adenocarcinoma cell line 

H23 cells were plated at a density of 5x 10 5 in four 
10-cm tissue culture dishes. Cells were transduced with 
50 MOl of the control Ad-CMV or Ad-CMV -RB2/pl 30 
and harvested after 48 h. Two tissue culture dishes were 
used to extract mRNA. One tissue culture dish was used 
to extract the proteins and one for FACS analysis. 

Northern blot analysis of samples transduced with 
RB2/pl30 showed an increased expression of the RB2I 
pl30 transcript more than 20-fold with respect to the 
control (Figure la). Western blot analysis showed more 
than 100-fold enhanced expression of pRb2/p!30 in the 
Ad-C M V- RB2\p 1 30 transduced cells (Figure lb). 
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Figure 1 Adcnovirus-mcdiated ovcrcxprcssion of RB2/pl30. 
Northern blot ami Wcsicrn bloi analyses (a. b) of RB2/pl30 in 
H23. H23-Ad-CMV. and H23 Ad-CMV-RB:/pl30 NSCLC cell 
line 
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Figure 2 Effects of pRb2/p 1 30-cnhanccd expression in H23 cells. 
FACS analvsis of H23-Ad-CMV (a) and H23-Ad-CMV-Rb2/pl30 
(b) infcctcdcclls. Rb2/pl30 ovcrcxprcssion resulted in 81.2% of the 
cells beinE accumulated in the G0/G1 phase of ihc cell cycle when 
compared" to the empty adenovirus (54%). The analysis was 
performed in triplicates with comparable results 



To confirm the effects of pRb2/pl30-enhanced 
expression in H23 cells, we performed FACS analysis. 
Fieure 2 shows that adenoviral Rb2/pl30 transduction 
resulted in 81.2% of the cells accumulated in the GO/ 
G 1 phase of the cell cycle when compared to the control 
(54%). z 

Oligonucleotide microarray assay following enhanced 
expression ofpRb2lpl30 in a human lung adenocarcinoma 
cell line 

H23 cells were transduced with Ad-CMV or Ad-CMV- 
RB2fpl30. After 48 h, 20 Mg of DNA-free total RNA 
from H23. H23-Ad-CMV, or H23-Ad-CMV-/?£2//>/J0 
cells was reverse transcribed and the double-strand 
cDNA was used as a template to generate Cy3-labeled 
cRNA probes and then hybridized co the Mergen 
H05000 oligonucleotide-based microarray containing 
more than 3200 genes that are involved in a variety of 
different processes. Analysis was performed by Mergen 
Ltd (San Leandro, CA, USA). Microarray experiments 
were performed comparing H23 vs H23-Ad-CMV. H23- 
Ad-CMV vs H23-Ad-CMV-/?*2//?m and H23 vs 
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H23-Ad-CMV-/v£2//>/jtf ceils. Duplicate experiments 
were carried out on a single total RNA preparation 
from the cells. 

In this study. 40 genes were downregulated more 
than 2.0-fold (Table 1). Figure 3 shows the plots of 
the differential expression of 3263 genes in H23- 
Ad-CMV vs H21-Ad~CMV-RB2lpU0 cells and H23 
rs H23-Ad-CMV-W2/p/J0 cells. Overall, the expres- 
sion of the majority of the spotted genes was not 
altered by RB2/pl30. Modulated genes were classified in 
Table 2 on the basis of a well-documented and 
established biological or pathological function of the 
encoded protein. The genes downregulated by pRb2/ 
pl30-enhanced expression belong mainly to the follow- 
ing categories: cell division; signaling and communica- 
tion: cell structure and motility; gene expression; 
metabolism; and disease. 



Validation of the oligonucleotide microarray assay using 
semiquantitative RT-PCR and Western blot analysts 
To determine the aene expression level of specific P Rb2/ 
pi 30 tareet eenes, semiquantitative RT-PCR analysis 
was used? A panel of 1 1 genes, randomly selected among 
the 40 identified bv microarray analysis, was analysed. 
We confirmed this by RT-PCR down regulation of * 
XfYB Cyc B2. Cvc Di t GRPR. KPNA2. MKK3. Nik. 
PCN4 PIM. PLK. and RAF- 1 (Figure 4). Genes highly 
downreaulated (range between 6- 17-fpld) in 

microarrav analysis such as PCNA. MKK3 B-tVfYB, 
and NIK showed a comparative downregulation in 
semiquantitative RT-PCR analysis between 7- and 30- 
fold. Genes still downregulated in microarray analysis, 
but at a lower extent such as RAF- 1. PIML CycDl 
GRPR KPNA2. and CvcBl* showed a comparable 



Table 1 Downrcgulaicd gene* by Rb2/pl3Q adcnovirus-cnhanccd expression 
GenBaiikTM ID Gene description 



Gene symbol Ratio I Ratio 2 Average 



1 


Ml 5796 




L36719 


3 


X13293 


4 


Z367I4 


5 


AF053306 


6 


LI9559 


7 


Y10256 


3 


D 14673 


9 


U33849 


10 


AF002322 


II 


LI 6862 


12 


L200IO 


13 


AB002359 


14 


X63692 


1 5 


U09559 


16 


AF011463 


17 


U2IS47 


18 


F046073 


19 


M25269 


20 


073335 


21 


X5I804 


22 


AB018330 


23 


L3380I 


24 


AFO0I9O3 


25 


M31735 


26 


M 32 110 


27 


X97630 


23 


M7348I 


29 


AF035536 


30 


D 10704 


31 


U 10339 


32 


Y 13936 


33 


AB000509 


34 


AJ0050I6 


35 


U6332 


36 


M54915 


37 


X59798 


33 


M63256 


39 


D26599 


40 


X03434 



Proliferating cell nuclear antigen 

Mitogcn-activaicd protein kinase kinase 3 

V-myb avian myeloblastosis viral oncogene homolog-likc 2 

c y clin F . x it 

Budding uninhibited by bcnzimidazolcs I (yeast homolog). ft 

Polo (Drosophia)-likc kinase 

Mitogcn-activaicd protein kinase kinase kinase kinase 4 
Kincsin-likc 2 

Proprotcin convcrtasc subtilisin/kcxin type 7 
Cyclin B2 

G protein-coupled receptor kinase 6 

Host cell factor CI 

FG A R amidotransfcrasc 

DNA (cytosinc-5-)-rncihyliransfcrasc 1 

Karyophcrin s 2 

Serine/threonine kinase 15 

TGFB-inuuciblc early growth response 

Budding uninhibited by benzimidazoics 1 (yeast homolog) 

ELK I. member of ETS oncogene family 

Uridine monophosphate kinase 

Putative receptor protein 

Calcium, calmodulin-dcpendcnt protein kinase kinase ft 
Glycogen synthase kinase 3 (I 

L-3-hvdro.xyacyl-Cocnzymc A dehydrogenase, short chain 
Polymerase (DNA directed). *> I. catalytic subunk 
Nucleolar protein 1 
ELKL motif kinase 

Gastrin-rclcusing peptide receptor , 
X-ray repair complementing defective repair in Chinese hamster ecus j 
Choline kinase 

Melanoma antigen, family A,3«'6 • 

Protein phosphatase IG (formerly 2C). magnesium-dependent, gamma .solorm 

TNF receptor- associated factor 5 

ATP-binding cassette, sub-family F (GCN20). member * 

TEA domain family member 4 

Pirn- 1 oncogene 

Cyclin Dl 

Cerebellar dcgcncration-rclatcd protein 
Protcasomc subunit. 0 type. 2 

V-raM murine leukemia viral oncogene homolog 1 ^ 



PCNA 

MKK.3 

B-MYB 

CCNF 

BUBIB 

PLK 

NIK 

KNSL2 

PCSK7 

CCNB2 

GPRK6 

HCFCl 

PFAS 

DNMT1 

SCPNA2 

STK15 

TIEG 

BUBI 

ELK1 

UMPK 

PMI 

CAMK.K.2 

GSK3B 

HADHSC 

POLDl 

NOLI 

EMKI 

GRP-R 

XRCC3 

CHK. 

VI AGE A 3/6 

PPMIG 

TRAF5 

ABCF2 

TEAD4 

PIMl 

CCS Dl 

CDR2 

PSMB2 

RAF1 



16.4 


17.5 


16.9 


12.4 


19.2 


15.S 


3.2 


13 


10.6 


6.6 


8.7 


7.6 


S.9 


6 


7.4 


6.5 


7.5 


7 


5.6 


7.6 


6.6 


5 


6.7 


5.3 


4.4 


6.5 


54 


5.3 


4,2 


4.7 


3.5 


5.7 


4.6 


3.2 


5.9 


4.5 


3.7 


5 


4.3 


3.5 


4.9 


4.2 


3.5 


4.1 


3.3 


3.7 


4 


3.3 


3.3 


4.3 


3.S 


4.3 


3.1 


37 


3.1 


3.1 


3.1 


3.2 


3.1 


3.1 


2.4 


3.7 


3 


2.1 


3.7 


2.9 


3.2 


2.6 


- 2.9 


3.2 


2.6 


2.9 


2.5 


3.3 


-- 2.9 


2.7 . 


3 


• 2.S 


2.4 


3.3 


2.3 


2.7 


2.7 


2.7 


2.5 


3 


2.7 


3 


2.3 


2.6 


2.1 


3 


2.5 


2.3 


2.3 


2.5 


. 2.2 


2.7 


2.4 


2.4 


2.5 


2.4 


2.5 




2.3 


2.2 


2.3 


2.2 


11 


2.1 


2.1 


2.2 


2.1 


2.1 


2.2 


2.1 


2.1 


2.3 




2.1 



Genes that are downregulated more than 2,.a>ld in^ponse to ^ 

identified as unique as mentioned in the Gen Bank -and « ^ ^^^j^^w^ fold of repression for each gene as 
determined by microarrav analysis of H23-AdCMV vs H23-Ad-CMV-RB2/pl30. Ratio - inaica o 
determined by microarray analysis ofH23 vs H23-Ad-CMV.RB2/pl30 
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Figure 3 Global comparison of gene expression in H23 vs H23* 
Ad-CM V and H23-Ad-CM V vs H23-Ad-CMV-RB2 cells. Each thi 
corresponds to the Cy3 fluorescent intensity of one single clement 
on the oligonucleotide microarray. A twofold change in expression 
is indicated with parallel lines marked as 2 



downresulation in semiquantitative RT-PCR analysis 
between 3.4- and 2.0-fold. PLK. which showed a high 
downregulation ratio in microarray analysis, failed to be 
validated by semiquantitative RT-PCR. In fact. PLK 
showed almost a twofold difference expression level by 
RT-PCR. 

Of the 1 1 transcriptionally downregulated genes that 
were studied by RT-PCR analysis, onlv seven genes (B- 
XfYB. KPNA2, MKK3. NIK. PLK, and RAF-I) were 
found expressed by Western blot analysis at a lower 
•level upon enhanced pRb2/p!30 expression with a ratio 
between 1.9- and 3.0-fold (Figure 5). As the MAGE gene 
family is composed of 23 related genes divided into four 
clusters and the MAGE- A subfamily comprises 12 genes 
highly identical in their coding sequence, we were not 
able to perform RT-PCR on this gene family, but we 
could confirm by Western blot analysis, the contingent 
downregulation of MAGEA-3/6 to enhanced pRb2/p!30 
expression. Surprisingly. PCNA that was highly down- 
regulated in the microarray analysis, also appearing 
modulated in RT-PCR, showed no protein expression 



Tahlc2 Classification of RB2/pl30-rcprcsscd genes by category 
~ ~~ Caws 



Catexftry 



ATPase/GTPasc/ATP binding/GTP binding 



abcf: 

KNSL2 

Calcium/poiassium/sodium/iron binding protein CAMKK2 



Cell cydc/cyclins 



Cell surface/antigen 



Chromosome/chroma tin/his tone 



Cytokines and growth factors 



Diffcrcniiation/dcvclopmcnt 



Diseases 



DNA binding/damage/ recombination 



G protein/regulators of G protein signaling 



Hydrolasc/hydrolysis/hydrolyscs 



BUBl 

BUBIB 

CCNBI 

CCNB2 

CCNDI 

CCNF 

B-MYB 

NOLI 

PCNA 

PLK. 

PPM1G 

ABCF2 

CDR2 

GPRK.6 

GRP-R 

KNSL2 

MAGE- A 3/6 

PCNA 

PMI 

DNMTl 

PLK 

XRCC3 

GRP-R 
PCSK7 
TIEG 
TRAF5 



Cvioskclcion/microiubulcs^micronlamcnis/moiiliiy CAMKK2 

EMK1 
KNSL2 



SUBl 

GSK3B 

B-MYB 

NOLI 

PIMI 

PLK 

RAFI 

TEAD4 

TIEG 

B-MYB 
CCNCH- 
CDR2 
ELK1 

MAGE-A 3/6 

NOL! 

PCSK7 

PIMI 

RAFI 

STK15 

DNMTl 

POLDI 

XRCC3 

GPRK.6 
GRP-R 
PMI 

ABCF2 
PCSK.7 
PPMIG 



Oncogene 
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Table 2 Continued 



Cutvgttry 



Kinases 



Lipoproteins/lipids 
Membrane trafficking 

Mitochondrial proteins 

Nuclear receptors/receptors 



Oncogenes 



Phosphatase/ protcascs/peptidase 



Signal transduction 



Synthetase/synthase 
Transcription/transcription factor 



Transporters 
Transferases 



Gaivs 



BUDI 
BUBIB 

CAMKK2 

CHK 

ELKI 

EMK.I 

GPKK.6 

GSK3B 

MK.K3 

NIK 

PLK 

RAFi 

STKI5 

UMPK 

CHK 

DNMTI 
KPNA2 

ABCF2 
HADHSC 

BUBI 

DNMTI 

ELKI 

GPRK6 

GRP-R 

NOLI 

PCNA 

PM1 

POLDI 

PPMIG 

T1EG 

TRAF5 

B-MYB 
ELKI 
P1MI 
RAFI 

PCSK7 
PPMIG 
PSMB2 

CAMKK2 

GPRK6 

GRP-R 

MKK3 

NIK 

PMI 

TRAF5 

GSK3B 
PFAS 

B-MYB 

CDR2 

ELKI 

HCFCl 

TEAD4 

TIEG 

ABCF2 
DNMTI 
PFAS 



Classifies the analysed genes on the basis of established biological or 
pathological functions of the eneoded proteins. Genes that arc listed in 
one category arc indicated in bold 
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HPRT 



Figure 4 Validation of oligonucleotide microarray results of 1 1 
selected genes by semiquantitative RT-PCR. RT-PCRs were 
performed using DNAsc- treated total RNA of H23. H23-Ad- 
CMV and H23-Ad-CMV.Rb2/pl30 NSCLC cell line. Amplified 
fniemcms of B-MYB (194 bp). Cyc B2 (217 bp), Cyc Dl (463 bp). 
GRPR (377 bp), K.PNA2 (304 bp). MKK3 (2l9bp). NIK (317 bp). 
PCNA (420 bp). PIMl (324bp). PLK (154bp). and RAFI (280 bp) 
sencs arc indicated. ACT-0 (626 bp) and HPRT (349 bp) genes 
were used as internal controls and were amplified from the same 
samples 



chanaes upon enhanced pRb2/pl30 expression. How- 
ever.lt has been shown that PCNA has a relatively long 
half-life that can extend beyond the S phase iruo the 
M phase and beyond into the GO phase of cells in rapidly 
Droliferating tumors. 



Discussion 

Durins the last 6 years, different studies demonstrated 
the involvement of RB2lpi30 in lung cancer. We have 
shown an inverse correlation between histological 
orading of lung cancer and the expression pattern of 
i*he tumor suppressor gene pRb2/pl30 (Baldi et aL 
1996, 1997: Minimo et aL 1999). Nevertheless the 
molecular mechanisms by which Rb2/pl30 participate 
in lung cancer are not yet clear. 

In the recent past, our group demonstrated that 
virally enhanced expression ofpRb2/p)30 caused about 



Oncogene 



# // / / 



B-MYB Z 



* 1 " ^ NIK 



E2F-1 
KPNA2 

MAGE-A 3/6 

MKK3 



PCNA 



PLK 



RAF1 



HSP70 



Figure 5 Validation of oligonucleotide microarray data by 
Western blot analysis. In all, 100/ijj of protein extracts from 
H23, H23-Ad-CMV. and H23-Ad-CMV-Rb2/pl30 cells were 
loaded onto SDS-PAGE eels and immunobloltcd with antibodies 
anti B-MYB, .E2F-1, MAGEA-3/6, MKK3. NIK. PCNA. PLK. 
and RAFl. Anti-HSP-70 was used as an internal control. The 
analysis was performed in duplicates with comparable results 



a 3.2-lbld or 69% tumor mass reduction in nude mice. 
These data strongly suggest the gene therapeutical 
potential of RB2/pl30 in lung cancer. In this view, the 
identification of genes modulated by enhanced pRb2/ 
pi 30 expression could serve as a valuable tool in 
identifying new targets for gene therapy in lung cancer 
as weil as to understand the molecular mechanisms 
behind this approach better. For these reasons, we 
decided to analyse the NSCLC H23 cells transduced 
with pRb2/pl30. 

In our study, we induced an overexpression on the 
RB2'tpl3Q gene in H23 cells (human lung adenocarci- 
noma) using an adenoviral system and performed an 
oligonucleotide microarray analysts to measure and 
evaluate the gene's expression profile. We identified 40 
genes that were downregulated more than 2.0-fold by 
pRb2/pl30 adeno vims-mediated overexpression. Regu- 
lated genes were summarized in 24 categories as shown 
in Table 2. We applied a cutoff ratio of 2.0 that has been 
commonly used in many microarray data analyses 
previously developed. We confirmed the modulation of 
11 selected genes by semiquantitative RT-PCR and ol 
six by Western blot analysis (Figures 4 and 5). 

B-MYB in our microarray analysis showed an almost 
10-fold reduction of expression with respect to the 
control 48 h after pRb2/pl30 transduction. Downregu- 
lation of B-MYB was confirmed by RT-PCR and 
Western blot analysis. B-MYB is a member of a family 
of transcription factors, which include c-MYB and A- 
MYB, that interact with the promoters of different genes 
involved in differentiation as well as cell cycle (Lyon 
et aL, 1994). It has been demonstrated that B-Myb can 



activate promoters, both via DNA binding-dependent 
and -independent mechanisms, following phosphoryla- 
tion mediated by cyclin A/E-cdk2 kinase (Lane et aL, 
1997; Saia et <//./l997). B-MYB is a positive regulatorof 
cell proliferation and its expression is required for the 
growth and survival of cells of different origins. B-MYB 
expression is stringently coupled to the proliferative 
state of the cell, and its overexpression promotes cell 
proliferation (Sala et aL. 1996). It has been reported that 
pRb2/pl30 increases during neuroblastoma differentia- 
tion with a parallel decrease of B-Myb (Raschella et aL, 
1997). Additionally, there is a large body of evidence 
that the negative regulation of B-MYB takes place, at 
least in part, at a transcriptional level (Raschella et aL, 
1996), through an E2F site in the promoter region of B- 
MYB when~RB family proteins form complexes with 
E2F factors bound to DNA (Zwicker et aL. 1996). 

Recently, it was also demonstrated that B-Myb is . 
abundantly expressed in primary lung cancers, suggest- 
ing that its overexpression may contribute to lung 
tumorigenesis (Hibi et aL, 1998). Our microarray data 
corroborating those of RT-PCR and Western blotting 
are in accordance with previously published studies. 
. Cyclin B2 (Cyc B2) in our microarray analysis showed 
a 4.7-fold reduction of expression with respect to the 
control 48 h after pRb2/pl30 transduction. Downregu- 
lation of Cyc B2 was confirmed by RT-PCR. Cyc B2 is 
a member ofthe cyclin family and a component of the B 
type. Recent data suggest that overexpression of human 
cyc B2, not due to mere gene amplification, could 
contribute through an alteration of the spindle check- 
point and chromosomal segregation, and to chromoso- 
mal instability of cancer cells, which is a phenotype 
persistently present in human cancers including lung 
cancer (Masuda and Takahashi. 2002). 

Cyclin Dl (Cyc Dl) showed a 2.1 -fold reduction of 
expression with respect to the control 48 h after pRb2/ 
pi 30 transduction by microarray analysis. Downreguia- 
tion of Cyc Dl was confirmed by RT-PCR. Cyc Dl is a 
cell-cycle regulator essential for Gl phase progression 
and a candidate proto-oncogene implicated in the 
pathogenesis of several human tumor types. The Cyc 
Dl gene maps to one ofthe most frequently amplified 
chromosomal regions (Hql3) in human carcinomas 
(Motokura and Arnold. 1993) and it is frequently 
amplified in breast and lung cancer (Keyomarsi and 
Pardee, 1993). Additionally, Cyc Dl is frequently 
overexpressed in the absence of genetic amplification 
in breast cancers (Gillett et aL, 1994) and its over- 
expression is associated with advanced local invasion 
and the presence of lymph node metastases in head and 
neck carcinomas (Muller et aL, 1997). Cyc Dl is 
amplified and overexpressed in NSCLC (Reissmann 
et aL. 1999) and its expression has been shown to be a 
negative prognostic marker in lung cancer (Caputi et aL. 
1999). confirming the hypothesis of Schauer et aL 
(1994). In fact, it has been suggested that lung tumors 
may evade cell-cycle control through abnormal expres- 
sion of cyc Dl. Additionally, cyc Dl overexpression and 
lack of expression of retinoblastoma protein is fre-- 
quently seen in lung cancer. 

Oncogene 



In our microarray. E2F-I was one of the genes that, 
unfortunately, resulted in not being statistically signifi- 
cant. However, because recently it was shown thai E2F- 
/ can act either as an oncogene or as a tumor suppressor 
gene and that E2F-1 overexpression may contribute to 
the development of NSCLCs by promoting proliferation 
(Gorgoulis et aL. 2002). we performed a promoter 
analysis of all the 1 1 selected genes using the TESS 
software searching for E2F-I and B-Myb DNA-binding 
regulatory regions. We found that the GRPR. MAGE- A 
3/6. MKK3. RAFL B-MYB, PLK. PCNA, PIML Cyc 
52 and Cyc Dl genes contain E2F-1 and/or B-Myb 
DNA-binding sites in their promoter, sequences. All 
these data strongly validate and explain the down- 
regulation of GRPR. MAGE- A 316. MKK3. RAF!. B- 
MYB. PLK, PCNA. PIML Cyc B2. and Cyc Dl in our 
system. In fact, the enhanced expression of pRb2/pl30 
downregulates E2F-1 and B-MYB expression (Figure 5). 
Interestingly, the promoter sequence of the NIK gene 
that does not contain E2F-1 and/or B-Myb DNA- 
binding sites is also downregulated by enhanced pRb2/ 
pi 30. No data are available for the human promoter 
sequence of the KPNA2 gene, instead. 

The gastrin-releasing factor receptor gene (GRP-R) 
showed a 2.7-fold reduction of expression with respect 
to the control 48 h after pRb2/p!30 transduction by 
microarray analysis. Downregulation of GRP-R was 
confirmed by RT-PCR. GRP-R is a member of the G 
protein-coupled receptor family and mediates important 
physiological actions of its specific ligand, the gastro- 
intestinaf hormone GRP. Many studies have found that 
GRP increases tumor cellular proliferation, leading to* 
the hypothesis that this peptide hormone is a mitogen 
important for the growth of various cancers (Rozengurt, 
1988). GRP has been implicated in the development of 
lung epithelium (Hoyt et aL 1993) and it is thought to 
be a growth factor involved in human lung carcinoma 
(Fa thief aL, 1996). GRP-R has mitogenic activity and it 
is produced in an autocrine fashion in small and non- 
small lung carcinoma (Siegfried et aL, 1999). Therefore, 
the preliminary data that pRb2/pl30 downregulates the 
transcription of this mitogen in lung cancer is an 
important hint that many and different pathways may 
be regulated by this member of the retinoblastoma 
familv. 

The karyopherin KPNA2 showed a 3.8-fold reduction 
of expression with respect to the control 48 h after 
pRb2/pl30 transduction by microarray analysis. Down- 
regulation of KPNA2 was confirmed by RT-PCR and 
Western blot analyses.. The active transport of proteins 
into the nucleus requires an array of proteins, including 
karyopherin. KPNA2 is one of six forms of karyopherin 
alpha proteins. Recently, it was shown that KPNA2 is 
able to interact with the p65 subunit of nuclear factor- 
k'B (NF-kB) (Cunningham et aL 2003). NF-kB activa- 
tion requires removal and degradation of its inhibitor 
kB, an event that occurs after phosphorylation of 
inhibitor fcB by a complex of inhibitor kB kinases. 
These events allow NF-k*B to translocate into the 
nucleus, to bind DNA, and regulate gene transcription*. 
NF-k-B nuclear translocation plays an important role' in 



preventing apoptotic cell death in some cancers. 
Recentlv,~ it was shown that inhibition of NF-kB 
enhances apoptosis in human lung adenocarcinoma 
cells in vitro (Milligan and Nopajaroonsri. 2001). 
Therefore, it is possible to hypothesize that down- 
regulation of KPNA2 could evade, in part. NF-kB 
nuclear translocation and subsequently, its antiapopto- 
tic effect. 

MAGE genes were initially identified because they 
encode tumor antigens that are recognized by cytolytic 
T lymphocvtes derived from blood lymphocytes of 
cancer patients (Van den Eynde and van der Bruggen. 
1997). The MAGE gene family is composed of 23 related 
genes divided into" four clusters. The MAGE-A sub- 
family is very appealing for its potential as antitumor 
immunotherapic because MAGE-A proteins are strictly 
tumor specific. Seven MAGE-A genes (MAGE-A L -A2, 
-A3, -A4 t -A6, -A 10, and -A 12) have been found highly 
transcribed in a large proportion of tumors of various 
histological origins, such as metastatic melanomas, 
testicular germ" cell tumors, lung cancer, sarcoma, 
mammary "tumors, and colon carcinomas (Ohman 
Forslund and Nordqvist, 2001). Recently, MAGE-A1, 
3. 6, 12, and 4b subfamily- members have also been 
proposed as candidate biomarkers in lung carcinoma 
(Sugita et aL, 2002). Our microarray analysis revealed 
thai MAGE-A 3/6 genes are downregulated with a ratio 
of 2.5 and these data were confirmed by Western blot 
analysis, suggesting an interesting link between MAGE- 
A 3/6 senes" RB2/pl30, and lung cancer. 

MKK3 showed a i 5.8-fold reduction of expression 
with respect to the control 48 h after pRb2/p!30 
transduction by microarray analysis. Downregulation 
of MKK3 was confirmed by RT-PCR and Western blot 
analyses. MK.K3 specifically activates p38-MAPK that 
regulates inflammation, apoptosis, and development 
(Schaeffer and Weber, 1999). Very recently, it has been 
shown that p38-MAPK is activated in NSCLC. 
suggesting an additional role of this kinase in malignant 
celTerowth or transformation (Greenberg et aL. 2002). 
Additionally. MKK3 has also been found to be 
homozygously deleted in a lung tumor cell line. All 
these data suggest that there is a close relation between 
this gene an cf lung carcinogenesis (Teng et aL. 2001). 

NIK showed a~6.6-fold reduction of expression with 
respect to the control 48 h after pRb2/pl30 transduction 
by microarray analysis. Downregulation of NIK was 
confirmed by RT-PCR and Western blot analyses. NIK 
is a protein kinase able to upregulate NF-kB (Malinin 
ei aL. 1997) that is highly expressed in NSCLC 
(Mukhopadhyay et aL, 1995), playing an important role 
in anchorage-independent and metastatic growth of 
lung carcinoma cells (Jiang et aL. 2001). 

The proto-oncogene Pirn- 1 {PfM-1) showed a 2.2-fold 
reduction of expression with respect to the control 48 h 
after pRb2/p!30 transduction by microarray analysis. 
Downregulation of PIM-l was confirmed by RT-PCR 
analysis." PIM-1 is a serine/threonine kinase that when 
overexpressed is involved in lymphomagenesis. PIM-1 
also associates with protein complexes necessary for 
mitosis and maps to an area on chromosome 6; which is 
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close to a region that is amplified in NSCLC (Bhutta- 
charya et aL 2002). 

Proliferating cell nuclear antigen (PCNA) showed a 
16 9-fold reduction of expression with respect to the 
control 48 h after pRblp!30 transduction by micro- 
array analysis. Downregulation of PCNA was validated 
only by RT-PCR. PCNA is a nuclear protein acting as a 
cofactor for DNA polymerase delta (Mathews et aL 
1984). PCNA has been used as a proliferation and 
prognostic marker in a variety of malignancies, includ- 
imf luns tumors (Kawai et aL 1994) and it is an 
important prognostic factor of clinical outcome in 
patients with NSCLC (Volm and Koomaei. 2000). 

Unfortunately. PCNA has not lived up to its early 
promise as a cellular proliferation marker. The relatively 
long half-life of the protein extends its expression 
beyond the S phase into the M phase and beyond again 
into the GO phase of cells in rapidly proliferating 
tumors. This could in part explain why we did not find . 
protein downregulation when compared to its transcrip- 
tional modulation following RB2/pl30-enhanced ex- 
pression. In fact, recently Markey et aL (2002) found 
that in a doxocycline-inducible RBI system. PCNA does 
not change at the protein level after 48 h from induced 
RBI expression, suggesting chat the failure to change 
could reflect a relatively long PCNA half-life. 

PLKl showed a sevenfold reduction of expression 
with respect to the control 48 h after pRb2/pl30 
transduction by microarray analysis. Downregulation 
of PLKl was confirmed by RT-PCR and Western blot 
analyses. PLKl is a member of PLKs serine-threonine 
kinases that are highly conserved during evolution. 
More evidence supports the concept that Plks regulate 
.different cell-cycle stages throughout mitosis, including 
its initiation by activating Cdc2 through Cdc25 and 
direct phosphorylation of cyclin Bl targeting Cdc2/ 
cyclin Bl to the nucleus. Moreover. Plks are key 
regulators of cytokinesis (Glover et aL, 1998: Nigg. 
1998). PLKl is involved in centrosome maturation. 
DNA damaae checkpoint adaptation, bipolar spindle 
formation, and activation of Cdcl6. Cdc27 as compo- 
nents of the anaphase-promoting complex for mitotic 
exit. PLKl is overexpressed in various human tumors, 
for example, NSCLC, squamous cell carcinomas of the 
head and neck, melanomas, endometrial, and ovarian 
carcinomas (Yuan et aL, 1997, 2002). It has also been 
recently shown that the constitutive expression of PLKl 
may contribute to cancer progression (Smith et aL. 
1997). PLKl overexpression is also a negative prognos- 
tic factor in patients suffering from NSCLC (Wolf et aL 
- 1997). All these data, taken together, suggest for the first 
time a link between a member of the PLK family. RB2/ 
p!30. and lung cancer. . 

RAF-l showed a 2.1-fold reduction of expression with 
respect to the control 48 h after pRb2/pl30 transduction 



by microarray analysis. Downregulation of WW was 
validated by RT-PCR and Western blot analyses. The 
serine/threonine kinase Raf-l is the first downstream 
effector of Ras to be identified. Activation of this kinase 
is necessary and sufficient for ^rt^***^ 
cascade consisting of MEKL2 and ^'"^ 
1994) Activation of the pathway RAFilMEK/tHK 
plavs an important role in gene expression control 
during cell cycle, apoptosis. cell differentiation and 
migration and is frequently associated with transforma- 
tion of primary cells and tumor progression. In act, 
recent studies suggest a possible role of Raf-l in lung 
tumorigenesis (Yano et aL 1999; Ramaknshna et aL 

2002). , f 

In a typical cell, it has been estimated that most ol the 
expressed genes can be expressed in abundance class at 
(or below) the detection threshold of microarray 
technology. This includes many genes that do not need 
to be expressed at extremely high levels in order to be 
biologically active, such as growth factors and tran- 
scriptional regulators (Fambrough et aL 1999). There- 
fore a more sensitive method such as semiquantitative 
RT-PCR analysis could be an acceptable strategy, not 
only to validate microarray data but also to identify . 
significantly induced genes whose transcripts belong to 
lower abundance classes. 

Our microarray analysis indicates . that adenovtrus- 
mediated enhanced expression of the tumor suppressor 
oene RB2/pl30 in the H23 human lung cancer cell line 
downregulated a variety of genes involved in many 
cellular processes including cell division, cell signaling/ 
cell communication, cell structure/mo tility, gene expres- 
sion, and metabolism. Adenovirus-mediated expression 
of pRb2/pl30 upreaulated a cluster of 28 genes with an 
average ratio between 2.0 and 8.5 (data not shown), and 
some of them have been found to be previously linked to 
lun° cancer. Further studies are needed to investigate 
the~$ignificance of these genes' regulation by pRb2/pl30 

better. , . nra 

This study identifies a cluster of genes that are 
modulated bv Rb2/ol30 expression. Although part of 
these senes could not be the target of pRb2/pl30 at its 
physiological level, most of them might mediate new 
potential*" therapeutical effects of RB2/pI30 in lung 
cancer. Therefore, we feel that these studies, by bringing 
about a better understanding of ample spectra ot 
proteins' expression, are essential because they could 
identify new cancer biomarkers and could facilitate 
matchin* the appropriate therapies to lung cancer, 
thereby maximizing therapeutic efficacy and minimizing 
toxicity. 
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